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Abstract 

The survival probability of a particle diffusing in the two dimensional 
domain x > near a "windy cliff" at x = is investigated. The 
particle dies upon reaching the edge of the cliff. In addition to dif- 
fusion, the particle is influenced by a steady "wind shear" with velocity 
v(x,y) = v sign(y) x, i. e., no average bias either toward or away from 
the cliff. For this semi-infinite system, the particle survival probability 
decays with time as t _1//4 , compared to t -1 / 2 in the absence of wind. 
Scaling descriptions are developed to elucidate this behavior, as well 
as the survival probability within a semi-infinite strip of finite width 
\y\ < w with particle absorption at x = 0. The behavior in the strip 
geometry can be described in terms of Taylor diffusion, an approach 
which accounts for the crossover to the i -1 / 4 decay when the width of 
the strip diverges. Supporting numerical simulations of our analytical 
results are presented. 
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I. Introduction 



Consider a particle which diffuses in the semi-infinite planar domain (x > 0, y) and is 
absorbed when x = is reached. The line x = can be viewed as a "cliff" and absorption 
at x = corresponds to the particle falling to its death. For this system, it is well known 
that the particle survival probability, S(t), decays in time as [1] 

S(t) ~ . (1) 



Here xq is the initial distance from the particle to the cliff and D is the diffusion coefficient. 
In this article, we are interested in understanding the time dependence of S(t) when the 
diffusing particle also experiences a "wind shear", defined as the velocity field v(x, y) = vx 
for y > and v(x, y) = —vx for y < (Fig. 1). Our primary result is that, although there 
is no average bias either toward or away from the cliff, the survival probability decays as 
t -1 / 4 , compared to the t -1 / 2 decay in the absence of the bias. This result is contrary to the 
naive intuition of a faster decay, as wind shear enhances longitudinal (x) diffusion, which, 
from Eq. (1), should reduce the survival probability. More generally, our interest is in un- 
derstanding the interplay between macroscopically heterogeneous convection and diffusion 
on first-passage phenomena. The wind shear geometry is a relatively simple example of 
such a system. In spite of this simplicity, relatively unusual first-passage characteristics 
occur. 

Our work is also complementary to a recent paper by Lee and Koplik [2] where the 
survival probability of a particle in the same wind shear is considered, but in a "two-layer" 
system with y unbounded, and with absorbing boundaries at x = and at x = L. For this 
system, Lee and Koplik gave asymptotic arguments to show that the survival probability 
decays as exp(— t/T) with T oc L. In comparison, for pure diffusion in the same domain, 
S(t) ~ exp(— t/T D ), but with a decay time T D oc L 2 . Thus, in accord with intuition, wind 
shear decreases the survival probability in the two-layer system (although the functional 
form of S(t) is unchanged). However, for the semi-infinite planar system, the wind shear 
has the opposite effect of enhancing the particle survival probability and moreover changes 
the decay exponent from 1/2 to 1/4. 



A fundamental aspect of diffusive motion in wind shear is that the particle spends 
relatively long periods of time exclusively in the region y > or y < before returning 
to y = 0. In fact, the distribution of time intervals between successive crossings of y = 
coincides with the first passage probability for a one-dimensional random walk to return 



long alternating flights where x increases or decreases linearly with time. For a particle 
which starts at the origin in an unbounded two-dimensional system, this structure for the 
longitudinal steps leads to the following unusual probability distribution in x at time t 



Rather surprisingly, the probability distribution is peaked as x — > ±vt and is a minimum 
for x = 0. These unusual features are one manifestation of the classical arcsine law for long 
leads in a one-dimensional random walk [3]. The interplay between these long excursions 
and the absorbing boundary are responsible for many of the intriguing features of the 
particle survival probability. 

In section II, we define microscopic lattice rules to model a particle moving in a wind 
shear. In section III, we present scaling arguments to determine the long-time behavior 
of S(t) for the semi-infinite planar system. Somewhat different behavior occurs if the 
underlying diffusion is only in the y (transverse) direction compared to the underlying 
diffusion being isotropic. While the former case is in some sense simpler, it does not have 
a smooth limit as bias vanishes. To understand this basic limit, we therefore introduce a 
variant of the original model which reduces to isotropic diffusion when the bias vanishes. In 
section IV, we apply the Taylor diffusion description to determine the survival probability 
in a finite width strip x > and \y\ < w, with reflection at \y\ = w and absorption 
at x = 0. Our results for the strip are combined with crossover arguments to provide 
additional insights into the corresponding first-passage properties of the planar system. A 
salient feature of the survival probability in the strip geometry is the non-trivial dependence 
of the survival probability on the initial distance from the particle to the cliff, the velocity, 
and the diffusion coefficient. We conclude and discuss some open questions in section V. 



to its starting point. Since this probability decays as t 3 / 2 , the particle has relatively 



[3,4]: 
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II. Wind Shear in Two Dimensions 



Our system is a square lattice with x > and y = denned to bisect the bonds 
which join the bottom row in the upper half plane to the top row in the lower half plane 
(Fig. 1). With this construction, there are no sites with y = 0, and ambiguities associated 
with assigning hopping rates from these sites are avoided. The hopping rates from any 
site in the upper (lower) half plane are spatially homogeneous. We consider two different 
hopping rules which correspond, respectively, to anisotropy and isotropy (for v — > 0) in 
the underlying diffusive motion. 

For anisotropic diffusion, to account for a variable transverse diffusivity and longitu- 
dinal velocity, we define the following hopping rates to six neighbors of a given site. For 
sites in the upper-half plane: 

Pl,±l = l^+^Dy, 

Pi,o = l(l + v)(l-D y ), 

P-i,±i = 4^ ~ v ) D y, 
P-i,o = \{l-v){l-D y ). 

Here the subscript on indicates vector displacement defined by the hop. In the lower 
half plane, similar hopping rates exist, except with an opposite sign for v. With these 
rules, there is a diffusivity in the re-direction whose magnitude depends on v and vanishes 
as t; — > 1. For the case of unit velocity and unit transverse diffusivity, these hopping 
rules reduce to a 2-site neighborhood. This appears to be the simplest implementation for 
random walk modeling of diffusion in a wind shear. 

While the above hopping rule is suitable for most of our purposes, it leads to patholog- 
ical behavior for v — > or for D y — > 0. Because useful insights can be gained by considering 
the crossover to purely diffusive behavior as v — > 0, we therefore define a second model in 
which the underlying diffusion is isotropic. For sites in the upper half plane, we define the 
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anisotropic diffusion (3a) 



hopping rates: 

Po,±i = D, 

v 

Pl,0 = D +~i 

y isotropic diffusion (36) 

P-1,0 = D- -, 

p ,o = 1 - 4D. 

Clearly, these rules can be applied only for D < 1/4 and v < 1/2 to ensure positivity of 
all the hopping rates. 

In the following section, we outline the asymptotic behavior of the survival probability 
with these two lattice realizations of wind shear. 

III. Survival Probability in the Planar System 

(a) Wind Shear with Anisotropic Diffusion 

Consider now the survival probability of a particle initially at (xo,0) in the semi- 
infinite planar system x > 0, with absorption at x = 0. From numerical exact enumeration 
of the probability distribution for the case v = D y = 1, it is evident that the survival 
probability decays as t -1 / 4 (Fig. 2). A relatively simple way to understand this result is to 
focus on those points where the particle trajectory crosses from y > to y < or vice versa. 
Since the transverse motion is a one-dimensional random walk, the probability distribution 
of times between successive crossings asymptotically varies as t~ 3 / 2 [1]. Consequently, the 
longitudinal displacement x(t) versus t is a Levy flight which consists of "segments" ti 
whose lengths are distributed according to the above distribution (Fig. 1) [5]. Here, the 
term "segment" refers to a connected portion of the trajectory with the y-coordinate having 
the same sign. 

The existence of a finite observation time t, however, implies that the segment length 
distribution is necessarily cut off at this time. With this cutoff, the average segment 
length is given by (t) oc /* t'-t'- 3 / 2 df oc t 1 ' 2 . Thus, as might be anticipated, a trajectory 
of t steps can typically be decomposed into Af oc \ft segments, each of length y/t. At 
the segment level, the probability that the walk does not reach x = is equivalent to 

5 



particle survival. Since there are Af independent segments, this no return probability 
should therefore vary as 1/VJ7 which, in turn, is proportional to t -1 / 4 , in accord with our 
observations. 

It is also instructive to consider the full dependence of S(t) on microscopic parameters. 
For our system, S(t) can only be a function of the basic variables x 0: v, D y , and t. Since 
S(t) is dimensionless, it is convenient to introduce the following groupings with units of 
time, T|| = xq/v and td = Xq/D v , and write for the survival probability 

*?(*)-/(-,—)■ (4) 

T|| T D 

A crucial fact is that the particle survival probability is governed by the difference in 
residence times within the regions y > and y < 0. As written in the original arcsine 
law (Eq. (2)), this difference is independent of D y . Combining this with the fact that 
S oc t -1 / 4 , the asymptotic form of S(t) reduces to 

*>«(?r-(sr- v 

To justify this rather unexpected behavior, it is instructive to consider the survival prob- 
ability on a finite width strip (Sec. IV). For a system with finite extent \y\ < w, then w 
and D y naturally appear in the combination rj_ = w 2 /D y . Since any finite value of D y 
leads to the same value of r± as w — > oo, it suggests that Eq. (5) should be independent 
of D y . Our numerical data supports this conclusion (Fig. 2), as long as T|| > 1/D y , so 
that there is mixing between the upper and lower half-planes, before significant absorp- 
tion occurs. However, our data typically exhibits a very small and unexplained residual 
systematic dependence of S(t) on D y . 

Another interesting aspect of the surviving particles is their spatial distribution (Fig. 3). 
Numerically, the surviving particles are predominantly within a region whose mean longi- 
tudinal and transverse positions are given by (x(t)) oc t and (y(t)) oc t 1 / 2 . These results 
can be justified heuristically. For walks to survive, there clearly must be a longer residence 
time in the upper half plane than the lower half plane. However, starting from (x,0), a 
walk can tolerate a short excursion into the lower half plane and still survive. If this excur- 
sion extends to a transverse distance y = — y/Dx/v, then the time required for the particle 
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to return to "safety" (y > 0) becomes of order of the time for the particle to be convected 
to the cliff. Thus if the particle enters the region y < —^jDx/v, it is likely to be absorbed 
(Fig. 1). Since surviving particles must almost always be in the region y > —^jDx/v, and 
hence predominantly in the upper half plane, this leads to the longitudinal displacement 
of the survivors being proportional to t. In a similar spirit, if the survivors are mostly in 
the upper half plane, then their transverse dispersion can be governed only by diffusion, 
so that (y(t)) oc t 1/2 . 

It is also worth noting that if (x(t)) oc t, then the typical value of y for which absorption 
occurs is —yjDt/v. Thus the two-dimensional system can be reduced to an effective one- 
dimensional transverse problem of a purely diffusing particle starting in the domain y > 
with an absorbing boundary at — yjDt/v. This latter system has been extensively studied 
[6]. It is known that S(t) oc t~ a with a dependent on D and v in such a way so that 
a ^ 1/2 for v — > and a — > for v — > oo. This connection between the wind shear and 
one-dimensional moving boundary problems provides further evidence of a decay exponent 
for S(t) in wind shear which is less that 1/2. 

(b) Wind Shear with Isotropic Diffusion 

When the underlying diffusion is isotropic, the survival probability again decays as 
t -1 / 4 in the long time limit (Fig. 4(a)). As might be anticipated, the effect of diffusion 
is subdominant with respect to convection in governing the value of exponent in the time 
dependence of S(t). However, for the isotropic system there is a crossover from diffusive 
behavior, for small v, to convective behavior, as v becomes large. By applying scaling to 
determine the nature of this crossover, we also determine the dependence of S(t) on the 
physical parameters of the system. This aspect merits emphasis, as there are important 
differences in the dependence of S (t) on system parameters with isotropic and anisotropic 
diffusion. 

For isotropic diffusion, the basic time scales of the system are 

2 

td = -jj, T|| = — , and t x = D/v 2 . (6a) 

These times are, respectively, the characteristic time to diffuse to the cliff, the time to con- 
vert to the cliff, and the crossover time beyond which convection dominates over diffusion. 
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It is convenient to measure these times in units of T\\ to give 

f D = Pe, f|| = 1, and f x = (66) 

where Pe = xqv/D is the Peclet number. In the case where Pe > 1, convection becomes 
established before the particle can reach the cliff and the behavior is the same as that 
discussed previously for anisotropic diffusion. On the other hand, for Pe < 1, the funda- 
mental times obey the inequalities td < f\ \ < f x . Consequently, particle death at early 
times is determined by diffusion and it is only for t > f x that convective behavior sets in 
(Fig. 4(b)). In dimensionless units, the early-time decay has the form S(t) ~ (Pe /f) 1 / 2 , 
while at late times the survival probability can be written as S(t — > oo) ~ A(Pe ) f -1 / 4 . 
Matching these two asymptotes at f x = 1/Pe fixes the amplitude A(Pe), from which 

S(t - oo) oc = D s/ H i/4 - ( 7 ) 

As shown in Fig. 4(a), the long-time data accords with scaling behavior predicted by 
Eq. (7). 

IV. Survival in the Semi-Infinite Strip 



It is instructive to examine the survival probability in the presence of wind shear in a 
semi-infinite strip of width \y\ < w. Because this system is effectively one dimensional, the 
longitudinal motion asymptotically reduces to diffusion in one dimension. However, this 
motion is properly described by Taylor diffusion which accounts for the interplay between 
diffusion and convection [7] . From this description, the form of the survival probability is, in 
principle, straightforward to deduce. It is then possible to infer properties of the survival 
probability in the planar system by letting w — > oo. However, there is an unexpected 
subtlety in the behavior of S(t) for the strip which depends on the relation between x 
and w. The resolution of this feature provides useful general insights about the nature of 
the survival probability. 

Taylor diffusion arises because the particle convects to the right while y > and 
convects to the left when y < 0. Since this switching between y > and y < is 
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governed by diffusion, the longitudinal motion is also diffusive. The characteristic time to 
switch between right-going and left-going segment is the same as the time needed to diffuse 
between the regions y > and y < (or vice versa), namely, tj_ ~ w 2 /D. Therefore the 
step length of the effective longitudinal random walk is i ~ v rj_ and the corresponding 
Taylor diffusivity is given by D\\ ~ £ 2 /t± ~ v 2 w 2 /D [7]. 

With this characterization of the longitudinal motion, we now investigate the survival 
probability in the strip in terms of the two basic time scales t\\ and tj_. First consider 
the limit t\\ 3> tj_, which can be re-expressed as xq ^> £. Physically, many longitudinal 
segments are needed before the cliff is reached, or equivalently, the walk encounters the sides 
of the strip many times before reaching x = 0. Therefore, before any absorption occurs, the 
trajectory has time to become truly one-dimensional. Thus we conclude that the particle 
survival probability is given by the suitably adapted one-dimensional expression, 

<?(t;x »£)~-^=. (8a) 

In terms of the basic time scales, this (dimensionless) survival probability can be re- 
expressed as 

S(t;x ^>l)~ -^=. (86) 

However, the converse limit t\\ <ri (xq <£.£), is more relevant for understanding the 
survival probability in the planar system with wind shear. In this case, there is significant 
loss of probability by absorption at x = before the sides of the strip play any role. 
Consequently, the survival probability should decay as (xo/vt) 1 / 4 at short times, as in the 
semi-infinite planar system, and cross over to a decay of the form Bt~ x l 2 when t 3> t± 
(Fig. 5). By matching these two asymptotes at t = rj_ we determine the amplitude B and 
thereby find for the survival probability, 

S(t;£ «^)oc (^) V4 - (9) 

A noteworthy feature of Eq. (8) and Eq. (9) is the opposite dependences in w. While 
S(t;xo 3> £) oc 1/w, S(t;xo <C £) oc w 1 / 2 (and is also proportional to x^ 4 ). In the latter 
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case, S(t; x$ <C £) must be increasing in w for there to be a slower time dependence in 
S(t; x <C £) when w = oo. The changeover from Eq. (8) to Eq. (9) as x becomes smaller 
than £ can be viewed as arising by xq "sticking" at £. This occurs because after one Taylor 
diffusion step, the walk is either absorbed or is reintroduced at a distance of order £ from 
the absorber. In the time t_i_ ~ w 2 /D required for this reintroduction, only a fraction 
(xo/wrjj 1 / 4 of the walks remain. Consequently, 

S(t; x < £) = S(t; x > £)\ Xo=i x 

Notice also that the limiting expressions for S(t) for x 3> £ and x G <ti £ both become of 
order w/y/Dt when a; = £■ This provides a useful self-consistency check on our approach. 

V. Discussion 




We have investigated the behavior of the survival probability for a diffusing particle 
with a planar absorbing boundary, or cliff, in the presence of a superimposed "wind shear" . 
Although this velocity field does not have any bias either toward or away from the cliff, the 
long-time behavior of the survival probability S(t) is strongly affected by the wind shear. 
Our primary result is that S(t) oc t" 1 / 4 for a semi-infinite plane, compared to a decay 
of t -1 / 2 for the survival probability in the absence of bias. Although our approaches are 
neither rigorous nor microscopic, they provide a good quantitative account of simulation 
results. It is worth noting that the image method, which yields the survival probability in 
the presence of a planar absorber for both diffusion and and uniform convection, does not 
appear to be generalizable to wind shear. 

Our prediction for the survival probability has applicability beyond the case of wind 
shear. Consider, for example, a particle diffusing with superimposed linear shear, v(x, y) = 
vyx, and with absorption at x = 0. For this system, exact enumeration of the probability 
distribution shows that the survival probability decays as t -1 / 4 [8]. This result is also sup- 
ported by the segment argument given in Sec. 111(a) which suggests that S(t) oc (rn/t) 1 / 4 , 
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where T||, now the typical time to convect from xq to the cliff in shear flow, is deter- 



mined by Vy/Dr\ \ t\\ = xq. Thus for the linear shear, we predict the following asymptotic 
dependence of the survival probability on model parameters, 




In fact, we believe that the t 1 / 4 decay should hold for any velocity field with v x {y) 



—v x (—y). For example, for power-law shear v(x, y) = vy \y\@ 1 x, the asymptotic form 



is expected. Notice, that only in the case of wind shear ((3 = 0) S(t) is independent of 
D. Another situation for which unusual behavior of S(t) can be anticipated is stratified 
random flow [9] in which v x (y) is a random zero-mean function of y. In this case, there 
could be different behavior for S(t) in a typical configuration of the velocity field and when 
averaged over all configurations of velocities. 

It would also be desirable to develop more rigorous and microscopic approaches to 
understand the first-passage characteristics of systems with various types of neutral bias 
fields. For the case of wind shear, the distinct nature of the problem for y > and y < 
suggests that the Wiener-Hopf technique may be suitable. In fact, this technique has been 
successfully applied to a related problem involving first passage in the presence of colored 
noise [10]. However, we have been unable to apply this method to our problem. 
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Figure Captions 

Figure 1. (a) Wind shear in two dimensions with a "cliff" at x = 0. The underlying lattice 
(thin lines) is defined with y = bisecting two horizontal rows of points. A typical 
particle trajectory in this system is sketched. If the particle strays below the curve 



y = —^/Dx/v, it is likely to be absorbed, (b) The same trajectory plotted with x 
as a function of t. The distribution of segment lengths U in this representation is 
proportional to t~ 3 / 2 (see text). 

Figure 2. Time dependence of S(t), based on exact enumeration of the spatial probability 
distribution for two "semi" -particles (each with weight 1/2) initially at (xo,yo) = 
(16, ±1/2), in a planar semi-infinite system with v = 1 and anisotropic diffusion with 
D y = 1, D y = 1/4, and D y = 1/16 (upper to lower data sets). The dashed straight 
line has slope —1/4. 

Figure 3. Time dependence of (x(t) n ) 1 / n and (y(t) n ) 1 / n for n = 1, 2 for the distribution of 
surviving walks for the case v = D y = 1 and xq = 1. The straight lines have slopes of 
1 and 1/2. 

Figure 4. (a) Representative results for the time dependence of the S(t) based on exact enumer- 
ation of the probability distribution for a planar semi-infinite system with isotropic 
diffusion and (i) {v,x ,D) = (0.2,1,0.25), or Pe = 0.8, (ii) (0.1,1,0.25), Pe = 0.4, 
(iii) (0.05,1,0.25) or Pe = 0.2 (lower to upper data sets, respectively). Plotted is 
S(t) Pe ~ 3 / 4 versus f. The dashed line has slope —1/4. (b) Schematic behavior for 
S(t) versus t on a double logarithmic scale with isotropic diffusion in the limit of 
Pe < 1. 

Figure 5. Schematic behavior for S(t) versus t on a double logarithmic scale for a finite strip of 
width w with anisotropic diffusion in the limits of (a) xq ^> i and (b) xq <^ i. 



13 



